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ABSTRACT 


The effect of cross-currents on the behavior of 
surface heated water discharges has been studied. Six 
experiments have been carried out with two runs in each of 
the large, moderate, and small Richardson number groups. 
Within the range of source Richardson numbers studied, the 
temperature profiles in both lateral and vertical directions 
have been found to be similar. Attempts have also been 
made to correlate the variations of axial excess temperature 
decay, jet half-width, half-depth, and trajectory with the 
ratio of the jet velocity to the velocity of the cross- 


current. 
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inlet area 

half-width of jet 

width of inlet 

average half-width 
upstream half-width 
downstream half-width 
temperature half-width 
velocity half-width 
half-width coefficient 
pressure drag coefficient 
friction drag coefficient 
half-depth coefficient 
temperature coefficient 
half-depth of jet 

initial jet depth 
velocity half-depth 
surface heat exchange coefficient 
pressure drag 
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discharge 

Reynold's number 

source Richardson number 
hydraulic radius 


temperature at any point 
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jet discharge temperature 
ambient temperature 

jet centerline temperature 
average excess centerline temperature 
jet velocity 

sur face jet velocity 

jet centerline velocity 

jet discharge velocity 
cross-current velocity 
lateral entrainment velocity 
Real entrainment velocity 
wind speed 

mean width 

jet thickness 

jet density 

ambient density 

wind stress 

entrainment coefficients 
velocity ratio 


Similarity parameters in vertical and 
horizontal direction respectively 


initial discharge angle between jet 
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CHAPTER« 4 


INTRODUCTION 


1.1 General 


The heated water discharge from the power generat- 
ing stations has certain effects on the environment. The 
effects are caused by the change in temperature of the lake 
or river which is the source of the condenser cooling water 
supply. 

The physical properties of water, such as density, 
viscosity, specific weight, surface tension, vapor pressure, 
gas solubility and diffusion, change with temperature. 
Among the above mentioned properties, the dissolved oxygen 
in water affects greatly the aquatic life. The excessive 
temperature increase in the water decreases the saturation 
levels of oxygen, and generally causes vertical stratifi- 
cation which cuts down oxygen circulation to the lower 
layers and also increases the evaporation. 

Biologically, many species of aquatic life are 
physiologically tolerant to a small change in temperature 
but will suffer greatly if temperature exceeds some 
Criticalilimit. “Loweistorms of aquatic lite, such as 
bacteria, algae, parasites and bacilli, can adapt to the 
change very easily and actually take advantage of it by 
multiplying rapidly. But some higher species of aquatic 


life, suffer most from the temperature change. Thus, the 
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population structure of the whole ecosystem in the water- 
body gets shifted or modified to an undesirable state. 

Heat addition to the waterbodies mainly comes 
from power plants, industry, municipal and irrigation 
return water. Parker and Krenkel (1969) indicate that 
‘power plants are responsible for about 75% of this waste 
heat discharge. The problem has been further aggravated 
by the increasing demand on electric power, the increasing 
size of the power plants, and the transition from fossil to 
nuclear fuel plants which have a comparatively lower 
thermal efficiency. 

1.2 Growth in electric energy requirement 

The consumption of electric energy has been 
rising sharply during the past decade as a result of a 
continuing development of the ‘electric society‘. Today, 
more household appliances are electrically powered. The 
growth in population, economy and the upgrading of living 
standards are the major factors contributing to this 
increasing demand. 

In Canada, per capita use of electric energy is 
the second highest among the nations of the world. The 
demand has tripled during the past 25 years. In 1972, the 
domestic demand for electricity has increased to 231,000 
million Kilowatt hours as compared to 40,000 million 
Kilowatt hours in 1945. In U.S., the situation is even 


worse as the demand is doubled every 10 years. 
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Figure 1.1 and Figure 1.2 show the forecast of 
electric energy demand in Canada (An Energy Policy of 
Canada, Phase..,- 1973) .and in U.S.A. (Parker and Krenkel, 
1969) up to the years 2000 and 2020 respectively. The 
projected figure indicates that 1,132,000 million Kilowatt 
hours of electric energy must be generated to meet the 
demand by the end of this century in Canada. 

Most of the power plants in Canada today are 
hydro or use fossil fuel. In 1970, hydro-power stations 
supplied 76% of total electricity generated in Canada and 
25% in U.S.A. Because sites suitable and economical to 
Revelop for hydro-power stations in North America have been 
used up and sites in remote areas are not favourable due 
to increased transmission costs, the demand for electric 
power would be met mainly by the fossil fuel-.and nuclear 
power plants in the future. The present forecast indicates 
that a decreasing proportion of electricity will be genera- 
ted from hydroelectric stations and an increasing propor- 
tion from nuclear plants by the end of this century. Hydro- 
electric energy production might double to about 310 billion 
kwh which is 30 % of the demand by the year 2000 but nuclear 
power plants will provide at least one-third of the total 
energy required. 

1.3 Cooling water requirements and efficiencies 

Power plants operate at a very low rate of thermal 

efficiency. A significant part of the heat is disposed in- 


to the waterways which causes thermal pollution. present 
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FIGURE 1.1 ELECTRIC ENERGY DEMAND PROJECTION IN CANADA UP TO 
THE YEAR 2000 
(AN ENERGY POLICY FOR CANADA, PHASE 1) 
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FIGURE 1.2 ELECTRIC ENERGY DEMAND PROJECTION IN U.S.A. UP TO 
THE YEAR 2020 
(PARKER & KRENKEL) 
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levels of efficiency are about 40% for fossil fuel plants 
and 30% for nuclear plants. The efficiency is controlled 
by thermodynamic considerations and by the economics of 
power plant design. Safety requirements and the limitation 
of material may not permit any substantial increase in 
efficiency and significant technological inne oerereriee on 
thermal efficiency are not foreseen in the near future. 

In the U.S., 50 trillion gallons of water were 
required for cooling purposes in 1968 and the amount is 
expected to be doubled by 1980. The report on Energy Policy 
of Canada gives the order of magnitude of known heated dis- 
Charges into Canadian water from various sources for the 
year 1970, which are summarised on table 1.1. It also 
gives the projected figures for the year 2000 based on the 
thermal efficiency mentioned. It is clear that if the 
current trends are to continue, the heated discharge, and 
the increased evaporation will have significant effects on 
our environment. 

1.4 Alternative Cooling Methods 

The waste heat from the power plants is usually 
released to the atmosphere or disposed into the waterways 
by one of the following cooling methods: 

(i) Once-through cooling system: This is so far the 
cheapest mode of heat disposal. Water is pumped 
from a lake or river and circulates through large 


heat exchangers. The warm water is then discharged 
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back into the source, most often by an open channel 
outfall. The excessive temperature is reduced by 
both direct heat transfer to the atmosphere and 
mixing of the warm water with cooler source water. 
There are two distinct modes of outfall discharges: 
a. Surface discharge from an effluent channel. 
b. Sub-surface discharge from submerged outfalls: 
1. Single jet and 
2. Multiple jet 
Because of the excessive temperature, the Warmer 
effluent is always buoyant relative to the source 
waterbody. Umrath (1971) favours the scheme of 
heated discharge as a surface jet. His argument 
is that if the jet discharged as a submerged jet, 
most of its heat is absorbed by the cooler water 
as it slowly rises to the surface thus affecting 
the aquatic life to a maximum. But if the heated 
water is discharged at the surface, most of the 
heat will be released to the atmosphere by evapor- 
ation and relatively little is transferred down to 
the cooler waterbody. In the latter mode of dis- 
charge, the impact on the environment and aquatic 
life by the heat discharge is reduced to a minimum. 
For economic reasons, once-through cooling systems 
are generally given first considerateion. 
Cooling Ponds: Cooling pond is one of the simplest 


methods of water cooling. Cooling of the warm 
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water occurs through evaporative cooling, direct 
transfer of heat to the air, and radiation. Low 
heat transfer rate and high land prices make this 
method unfavourable in some areas. Generally, one 
to two acres of land surface is required per 
megawatt of installed capacity. Berman (1971) 
suggests) a iheativtransfer rate ‘of7/5 vto ob50 
BTU/hr /£t?; thus, for power plants with heat load 
greater than 800,000 BTU/hr, acres of land must be 
acquired to achieve the cooling purpose. Besides 
the large area acquired, other disadvantages in- 
clude the requirement of blowdown of the concen- 
trated dissolved solids and the collection of 
impurities due to large surface open to the 
atmosphere. Local meteorologic conditions such 
aS air temperature, relative humidity, wind speed 
and solar radiation must be studied before the 
determination of the location and size of the pond. 
(iii) Spray pond: The addition of spray modules to the 
cooling pond can increase the heat transfer rate 
but also increases the water loss due to evapora- 
tion. The heat transfer rate is about 4,000 to 
10,000 BTU/hr /£t? for spray pond as compared to 
75-150 BTU/hr/ft* for cooling ponds without spray 
modules. Spray ponds, however, are usually un- 
suitable for a cooling range greater than 18°F. 


In general, it is also uneconomical for large size 
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plants. As the size of the pond increases, the air 
passing by the sprays is already saturated so the 
sprays modules on the downwind side are not very 
effective. 

Dry cooling towers: The basic idea of the dry 
tower cooling system derives from the cooling 
system of an automobile. The heated effluent is 
circulated through the radiator and the air is 
forced to pass the radiator by the fan. The 

water is cooled by heat transfer through the walls 
of the radiator to the air. This system can be used 
where fluids to be cooled are at a high temperature. 
It eliminates water problems such as availability, 
chemical treatment and corrosion; but the instal- 
lation and maintenance costs are high. 

Wet cooling towers: This system works like the 
spray pond but the heated effluent is sprayed 
inside a large tower. Major mode of heat dissi- 
pation is evaporation which can be achieved by 
movement of air Beck the water droplets. Air 
movement can be obtained by using fans or by 
natural convection. The cooled water collects in 
the basin at the bottom of the tower. There are 
many types of wet cooling towers. The type choosen 
depends upon the requirements of the particular 


system. 
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In general, the choice of cooling method depends 
on location, available technology, economy and the objectives 
PonBeraecained.” “tment cases, the once-through cooling 
system is preferred if the plant location is close to a 
water source because of economic reasons. Cooling ponds 
require too large an area and the spray ponds are not very 
iefrecuuventon largevolanus.) VeLy high capiual wand mainte— 
nance costs associated with cooling towers make this 
system unfavourable. Besides, localized ground fog 
problems associated with cooling ponds and cooling towers 
make the once-through cooling method preferable where 
feasible. 

Wao DESCription Of ‘the! problem _ 

The impact on environment by waste heat disposal 
has aroused unprecedented public awareness recently. Some 
states in U.S. have set. up regulations restricting the rise 
in temperature in the receiving water to some distance away 
from the point of discharge in order to protect the envi- 
ronment; but some states have no restrictions at all. [In 
Canada, investigations have been carried out in varies lakes 
to study the waste heat distribution and its effects on the 
ecological system. 

The discharge of heated water on the surface of 
a lake spreads out like a surface buoyant jet. The temper- 
ature distribution in the discharge is considered as function 


of the discharge channel geometry, the initial velocity and 
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temperature rise, the magnitude of the ambient current, 

and the surface heat transfer rate. The spread of the 

warm water has two distinct regions given as: 

(i) Near-field region or jet region where the momentum of 
heat discharge must be considered. 

(ii) The far region where surface heat transfer aspects 
and ambient turbulence of the receiving water body are 
Significant. 

From the design point of view, the jet regime is of 
particular interest because most of the criteria are set 
£Or within this region. The outlet must be so designed 
that it meets the prescribed temperature criteria. Also, 
a complete knowledge of the temperature distribution is 
important for the selection of the location of the intake 


in order to prevent recirculation problems. 
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CHAPTER 2 


LITERATURE REVIEW 


2.1 General 


The purpose of this study is to study the 
characteristics of heated water discharges into a large 
water source with small ambient current such as lakes. As 
the heated water is discharged into the ambient water, it 
creates a tangential separation surface due to the velocity 
difference between the jet and the ambient. This shear 
between the moving jet and the ambient fluid creates tur- 
bulence which travels both inwards to the jet centerline and 
outwards towards the ambient fluid. Eventually, this tur- 
bulent mixing action reaches the jet axis at some distance 
from the outfall and the jet centerline velocity begins to 
decrease. The region before the mixing action intrudes into 
the centerline is called the zone of flow establishment and 
is of rather limited extent. The region further downstream 
(in the near field) is called the zone of established flow 
and is generally assumed that the velocity profiles are fully 
developed and are similar. These two regions combined to- 
gether is known as the near-field region. At some distance 
from outfall, the kinetic energy of the discharge would be 
sufficiently dissipated; the ambient turbulence dominates 
over the jet turbulence and governs the behavior of the jet. 


This region is classified as the far-field region and is 
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generally much larger than the near-field region. In the 
presence of an ambient current, the jet is forced to bend 
over in the direction of the ambient current due to the 
pressure difference created between the upcurrent and lee- 
ward sides of the jet. The curvature of the bent-over jet 
may be large or small depending on the velocity ratio, 
buoyancy and other factors. The addition of cross-flow 
momentum due to entrainment is the other basic mechanism 
causing the plume to level off. As soon as the jet bends 
over and becomes aligned with the ambient current, the 
cross-current interaction ceases for all practical pur- 
poses. 

A review on the existing literature reveals that 
few investigations have been done on simulating the char- 
acteristics of heated water discharges into a large water 
source. Most investigators study a two-dimensional sur- 
face jet assuming that the entrainment in the vertical 
foe en is negligible. This assumption may not be jus- 
tified if the local Richardson number is in the small 
Richardson number group (i.e. less than 0.1) in which case 
the jet thickness increases in the longitudinal direction. 
Therefore, a study on the three-dimensional surface dis- 
Charge with ambient current is of immediate interest and 
the author hopes to provide some useful experimental in- 
formation towards solving this complicate problem in this 


study. 
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2.2 Integral Prediction Methods 


This method is the most common approach in the 
model development. A set of differential equations for 
conservation of heat, momentum, and volume are derived 
for a control volume, and then integrated over the jet in 
the vertical and lateral directions to obtain the general 
expressions. 

To solve these differential equations numeri- 
cally, some assumptions must be made. In addition to the 
Boussinesg Approximation, which states that the density 
difference can be neglected except in terms involving 
gravity, velocity and temperature profiles are assumed to 
be similar in consecutive transverse sections. Morton 
(1961) states that the assumption of similar profiles 
suppresses all details of the transverse structure of the 
jet, and any profile shape can be used without loss of 
additional physical information. Most researchers propose 
a Gaussian type distribution which generally agrees well 
with the experimental data. 

The entrainment velocity is usually related to 
some characteristic velocity of the jet by coefficients. 
Platten and Keffer (1968) proposed an expression for en- 
trainment velocity for nonbuoyant round jets in air as: 


(U_, - U. cos @) + ao U (cos © - cos 8.) (23224) 


i ls 


14 


ae, Ri ae) rit 
tes a me ; ae ry 
cA ea ey 7 
; rh om Ny 1 an 
PARE waaay tt 
v, : 
; Welt 
i mid } 
‘i : 
bewi xob a6. ! 


ol Jab edt “30 


. io ih a4 yy 


lexvenop orld nis3ce: oF protioestt ls2et fi de feoks ae ott 
j Al * 4 ki ge as i 
Fs | iy | 


ieee snokseupe isL¢nero2? ib am | 
abe od Anois Lbs ar eben ed tenn 

wstiensb EYE) dad? wegeta foidw cn ea A ponte ue 

pnivfovni amzet mi) tae oie fyadobiliead od neo ebae’ 33 mi 


od bemogse eve nefiiow oiids:sdmes bis vitbotey \yivedy | i 


ii —/ 


| 2 tous anayd eff? 20 el ~ 
eit 20 @ 4h CUTIE, SR toOvenats a @, Ra apa 
to scot Juansiw ete et har aqats, obi 


J 


nod 1eM ‘Pa hIHeE ga Taverner ovisuoes 


tn 


4 | 
ao fhio% raLiaba > eo naidqgmimas edt F 


Bodo ny ‘eieriotaoney Slab ottemzolal I 
ilew soo 7p ikeseoge tobi rot 


ot bedslex : esti hi 


ea Sanna 
“Songhai ees cd a ons 


ne “a De 


where 2, is the angle between direction of cross flow ve- 
locity and centerline of inlet; os and a are entrain- 
ment coefficients which have to be determined empirically. 
Ene ssecond term jin Equation 2.2.1 is included! to account 
for the entraining effect of internal rotation imparted 
to the jet by the cross-flowing current. A simpler form 
is adopted by Morton (1956), and Motz and Benedict (1972). 
They assume that the entrainment velocity is to be propor- 
tional to the magnitude of difference between some char- 
acteristic velocity in the jet and the velocity of ambient 


unaffected by the jet. This gives: 


cos @) (252.2) 
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WS ay (U,, aia tb 
which has exactly the same form as the first term in 
Equatwony 2.24 1.. 

To develop the governing equations, we consider 
the regimes of moderate and large Richardson number (i.e. 
with the source Richardson number greater than 0.1). In 
these two regimes, the thickness of the surface jet can be 
assumed to be approximately constant at a very short dis- 
tance from the outlet. This assumption constitutes a 
great simplication in analytical modeling. 
Referring to Figure 2.1 and 2.2, the derived 


equations for conservation of volume, momentum, and heat 


flux for a control volume can be written as: 
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Conservation of volume flux 


The rate of volume flux change in the jet be- 
tween two sections along the axis should equal to the 
ambient flow entrained within that two sections. This gives: 


d 


ee 
eV. (@) mee f (221253 ) 
ax, / 4 He U dy, dz = 2 ie Zz 


Oo 


where |'V is the depth-average entrainment velocity; by is 


el 
the distance measured from the axis towarding to the edge 
to the point where the velocity is 1% of the centerline 


velocity and is usually taken as 2.25 times the half-width 


(by) of the same section and Zo is “the depth of the jet. 


Assuming: 
Us aT 
aes £145) Set (2.5 Aa) 
U <1 
Bate = ce EET a f£ 
u. Ae pr) (2.2. 4b) 
gives U_ = 7 


where b,, and h,, are the half-width and the half-depth in the 
lateral and vertical directions respectively. Substituting 
Equation 2.2.5 into Equation 2.2.3 gives the final’ form of 


the conservation of volume flux as: 
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Conservation of Momentum Flux in X Direction 


The momentum flux of the jet in the X-direction is 
balanced by the pressure drag created by the cross flow 


in the same direction. This gives: 


Z on AAS Uno. ,(U, coss0@). dy... dz 
dx, 0 0 1 i 
2 2 
ran oe cos 6 . 
A ——+5;-—— Zz sin a C2257) 


Applying the relationship as in Equation 2.2.5 to Equation 


Zecisf leads to: 


2 2 
= da 2 WU ocosi ne 
2 dx, Ze b. ne COS, 0 =. hs) Zo Sin 6 (2.2.8) 
225 yO, 2s 5 5 
where F_ = I iP Cee) (£,( 6 )) duvde 
= constant 


note that the friction drag is neglected in the above 


equation. 
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Conservation of Momentum Flux _in Y direction 


The variation of the momentum flux in the Y 
direction can be written as: 


b Zz 


2 
feo. a a. ne 
fae ; sin © fe * 2 p|V.| ZS UL 


Ca Seer a oa Zz cos ® (2.2.9) 


The second term in Equation 2.2.9 is the momentum flux 
increased or decreased due to the entrained ambient fluid. 
Combining with Equation 2.2.5, the equation for the con- 


servation of momentum. flux in Y direction is: 


2 3 
2 U, cos @ 
: eo = ik 
u_ b hh sine 2 Vig Cae Ue icecageaea 


C2 De) 
wherein the direct contribution of wind shear stress on 


the water surface has been neglected. 
Conservation of Heat Flux 


Neglecting the heat loss on the surface, the 


equation for conservation of heat flux is: 
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where br is the location at which the temperature excess 
is about 1% of the centerline temperature excess of the 


Same section. Further simplification leads to: 


age RO Doh, S570 (2.2.12) 
aa Re UO 

where 

ATe= Ee - (2, 

Geometrical Equations 

Se = cos @ (2.2. 3a) 
A 

cS ae 

Ava ano 20 (22261 3D) 
1 


We now have four integral equations (Equation 2.2.6, 

De 25078207107 292.12) and@apalr® of ‘geometrical’ equations 
(Equation 2.2.13a and Equation 2.2.13b). To proceed with 
the solution of these equations, we have to make further 
assumptions regarding the velocity and temperature dis- 


tributions, and entrainment velocity. 


21 


t " i. :* 
we eds | 
i i 


(et «& am) 


: | Bes ee : 
ec “teat 


ar ene Mie“. 


ise? 4%. S) 


(d€L.8-.8) 


Velocity distribution 


UP Ue sins As ae 
—— expe onees (2+) 5) (262. 14a) 
a Uy Sin 6 b, 


Temperature distribution 


T - De Yj? 
ee exp’ (“0.693 "(——)*) 

b Deo ela 
Ter T ( ) 
Entrainment velocity 
VS a COR U, sin 8) (222445) 


Numerical methods can be used to solve the above differen- 
tial equations to yield the jet trajectory, velocity and 


temperature decay, jet half-width and half-depth. 
2.2.1 - Carter 4s) Model: 


Carter developes a mathematical model based on 
his experiments and the experimental results of Rouse on 
air jets. He treats the problem as a two-dimensional 
surface discharge with uniform cross-current and assumes 
top-hat distributions for velocity and temperature across 
the width and depth of the jet. He considers a jet ina 
cross-current as an obstruction in a uniform stream with 


some of the ambient current intruding into the jet and the 
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rest of the oncoming stream being deflected around the 
jet. 

According to Rouse and Carter, the faces of 
the inlet exert a force downstream due to the pressure 
gradient across them, causing the jet to attain momen- 
tum at the orifice. The momentum equation in the Y 
direction can be written as: 


az 
dent. cy bh.) sine = U, V, 0. ho dx, 


Z 
= kp, C Uy ho cos © dx, (226271. 1) 


d 
Balancing the rate of change of X momentum with the X 
component of the pressure drag yields the momentum 
equation in the X direction. 
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2 
aor le b) cos © = -4 a Cc, U, hy sin e dx, 


(2 oe leee) 


These two momentum equations are then integrated from 
x, = 0 to the arbitrary centerline position x) and then 
ainvading Equation 12.12) bw byrEquatiron.<2)..2\1 sbigives: tthe 


slope of the centerline trajectory at the point (x,y): 
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where f is defined as the fraction of the control volume 
that is ambient fluid. Carter derived the equation for 
jet trajectory for both the zone of flow establishment 
and zone of established flow. 

For the zone of flow establishment, Carter 
assumes no dilution, i.e. f£f = 0, which does not appear 
to be a correct assumption. With this assumption, 
Bouatdon 2.2.1.3 -1stsolved’to yield a circular jet 


trajectory described by: 


pra o (2021.4) 


For the region of established flow, Carter 
assumes U = ae This assumption basically replaces the 
need for a continuity equation but does not seem logical. 


Besides, there is no hard evidence to suppport this as- 


sumption other than the author's observation that the flow 
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turned rapidly paralled to the ambient current. Thus, 
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FOr l= f= constant, Carter States that. tne 
trajectory is an ellipse with the same center but with 
semi-major and semi-minor axis given by aa b /Cg and 
(1 - 6)? 2 ee DE) Ge respectively. These axes are proved 
to be incorrect later by Policastro and Tokar. The de- 
crease in temperature along the jet axis is related to 
dilution which Carter defines as the fraction of a unit 
volume that is initially part of the jet discharge; the 
remainder is part of the ambient stream and denoted by 


f. Their relationship can be expressed as: 


Tn - T AT 
i Ns sy hae pa aah ali i (2e2. 0.0) 
To - T, AT) 


Carter plots his experimental data on a log-log graph 
paper and shows that the equation 

x -i 
AT i" 1 4 


=|— (2e0218-7) 
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fits the data well. This equation also reveals that the 


excess centerline temperature decay has a form apparently 
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independent of the velocity ratio a. Data analysis 
gives a value 2.5 bo as the length of the temperature 
flow establishment region but 3 b, is believed to be more 
adequate by other researchers. 

The determination of 8, and Cay dependwheavily 
on the measurements made by Rouse on an air jet in an air 
cross-current. The aurher derives the expressron for 6, 
by balancing the pressure gradient in the Y direction with 
the rate of change of Y momentum across the orifice. The 


expression gives: 


jeat C 
SinsO. =o Wd Cia) 
oO O 2) 


The value of C, has been determined as a function ofa 
from Rouse's experimental data on boundary pressure dis- 
tribution. An independent estimate of Cg can be made 


from the measured length of the zone of flow establishment. 


Since the trajectory in the region of flow establishment is 


described by a circle, the arc of the circle is then equal 
to the central angle measured in radians multiplied by the 


radius. Thus: 


2 2 
2, oO bo | 27 | 
Zao De mI a (8. se eo.) 360 | (Oe 2a eo) 


ONg 


26 


(a. fOy8) 
6 he noteaah 28 henshiesab need raat oti 
<gih, Sandaeneh Seaibinkive 9, 6888 (donoae xo 

shan ott TBO -: Ro! Sasntdeo ie pre 
_aheinte i Lches a= wOEd ihe! was bits 20 Hehe, “ Wes mds ao 7 
ai drenmteitceves age to nokge: ed+ ok yagag ! 3 
Leupe mond e4 esovke ait, 3 aes) Or) ae 
ode Yet ee hina fomunnen,; { 


(@.t.%.8) 
; oa “a ui Se 
A | i f vy 
1 , ae “ : 4 
- y af i 7 ny - 
re Ge \ es os) aan of ual 


C22 els LO) 


where 8. is the angle at the end of the flow establishment 
region. Based on previous assumption that the Y velocity 
component in the region of established flow equals Uy 


this gives: 
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SOlMing ROUat LON 26 2aknOe (eee. LO, 2.2211, vields the 
values for a. and Cq- 
Carter's model is different from other models 

in several aspects. The drag force defined as: 

pecs A v7 

sa a 
Dp cal a 2th 2) 
p 2 
is an artificiality to match theoretical jet trajectories 
with experimental data. Fan, and Motz and Benedict adopt 


the same form but with an additional term. Fan chooses 
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whereas Motz and Benedict choose 
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so to cause Dy to approach zero as the jet orients itself 
parallel to the current. The length of flow establishment 
is assumed constant with respect to velocity ratio but 
this assumption is not supported with experimental data. 
The continuity equation is not used. The assumption that 
the Y component of jet velocity equals ambient velocity 
in established flow region is also a doubtful approxima- 


tion especially in the vicinity of the end of the flow- 


establishment region for large velocity ratios. 


2.2.2 Hoope's Model 


Hoopes et al (1971) derive a mathematical model 
including direct wind stress by applying mass, momentum, 
and heat conservation to a control volume of the jet. The 
model is two-dimensional and the Richardson numbers are 
less than 0.1 which belongs to the small Richardson number 
group in our, system ofjolassitication. 

To develop the conservation of mass equation, 
the flow across a surface of the control volume perpen- 
dicular to the jet axis is considered. The increase in 


jet flow due to entrainment is: 
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£(y,/b) is an assumed form for the lateral velocity dis- 
tribution. The authors also assume that transverse velo- 


city profiles are similar at all sections. This leads to 
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If n > b, then the evaluation of I, will give a value 


1 
1.49, s“Bquation, 2.2.2.3h4 and 2.2:2.3B describe = Gaussian 
distribution for y, < n and a constant velocity for ae 
The entrainment velocity is simply assumed to be propor- 


tional to the jet center-line velocity and the amount of 


fluid entrained into the control volume is: 
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- o Ax,) a, U, (2 25 ax) (2.2.2.4) 


Equating BQuation 2.2025 band 2.2.2.4 yields the continuity 


equation. 

5 heal W eats 8, 2 Oy 
m = U (2298995) 
Ax, 1a m 


where ly is the entrainment coefficient. The momentum 
flux 1s expressed in KX and Y direction separately.” In 


X-direction. 
EF, =F + F = ay (2b dx, ) Sin 8 O22) 


where Spe and ae are the net pressure and shear forces in 
X direction respectively. g is the angle between wind 


direction and Y axis. The wind stress may be expressed as 
T = p C W (22627) 


where W is the wind velocity. The authors consider only 
the wind shear at the water surface; the pressure force, 


drag force, and the centrifugal force are all neglected. 


Therefore 
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Teut(y,/b) is given by Equation 2.2.2.3, ‘then I, will have 
om tlie Or I. 20. ; SUubstleucing Equation! 2.222.7 and 
2.2.2.8 into 2.2.2.6 gives the momentum equation in the 


X-direction. 
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final form of the momentum equation in the Y-direction as: 
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Before solving the governing equations, the authors 


introduce three dimensionless variables. They are: 
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Wochevhe ard of Equation 222.2. 13A -to.2..2.2 13H, Equation 
Pree eee 2. 2 -anGe. 22. LZ Can. be solved) numerical ly 
to yield the characteristics of the jet, i.e. Une bye; 
in terms of dimensionless volume and momentum flux. These 


Characteristics equations are: 
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Applying conservation of thermal energy to the 


control volume gives: 
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in which K is the surface heat exchange coefficient. The 
authors assume the similarity profile for the lateral 


temperature distribution as: 
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if the surface heat loss is negligible, 1.6. K = 0, ther 


Equation 2.2.2.20 reduces to: 


ee ae ee is Ty (202) 
A weak point in this model involves the author's 
analysis of entrainment. Ellison and Turner, Keffer ond 
Baines, Monton, and Fan have all defined entrainment in 
terms of the difference between a characteristic jet 
velocity and the velocity of the moving ambient current 
instead of directly related to the centerline veloity 
as in Hoope's model. The assumption of similar profiles 
for velocity and temperature makes the model capable of 


handling only very small cross flows. They assume a linear 


spreading of the jet and this may not be true if the level 


of the ambient turbulence is low and buoyancy effects are 
predominant. The model also fails to handle the lack of 
Symmetry on windward and leeward sides of the jet which is 


reflected in an obvious change in entrainment coefficient. 


2.2.3 Motz and Benedict's, Model. 


This steady, turbulent, two-dimensional jet 
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model is developed based upon the framework established by 
Morton, Fan, Hoopes and Carter. The argument on no ver- 
tical entrainment is discussed by Ellison and Turner in 
Which they find that the Richardson number increases 
gradually to one or greater along the jet axis regardless 
the source Richardson number. When the Richardson number 
is close to or greater than unity, buoyancy forces domi- 
nate over inertia forces and entrainment of the ambient 
water occurs only at the sides of the jet and not to any 
Significant extent at the bottom. 

A control volume is segregated from the jet to 
develop the governing equations. The equation for the 
volume flux can be written by balancing the rate of 


Change of volume flux and the fluid entrained. 


Thus: 
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The entrainment velocity is defined as: 
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which is an improvement over the form used in the model 
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of Hoopes et al. The entrainment coefficient ay. as in 
other models, has to be determined empirically. The 
momentum equations in the X and Y directions are: 
JaENeGhon- 
Smeets cds.e)) r= eaggee COE Se epee 
dx, / T S) 
d oe 2V2 
oe es = % 
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i clei ree | 
ee (282 252) 
respectively. The heat conservation equation is simply 
expressed as: 
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and also the two geometric equations: 
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Thus, there are six equations which are functions of the 


jet axis distance x Before solving the above differen- 


L 
tial equations, the temperature and velocity ere nies 
have to be assumed and some dimensionless parameters are 
introduced. Motz and Benedict assume similarity in both 
the temperature and velocity profiles dnd’ choose an 


exponential form for the similarity profiles. The Gaussian 


profiles chosen are: 
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The dimensionless parameters are: 
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The drag force D_, is assumed to be related to the ambient 


velocity in the following form: 


a —. sin © (2at2e 55 L4) 


With the additional information on the drag and 
entrainment coefficient, the authors are able to solve 
BOMACLOM U2. 265.5 Lele 2526564) Ulaedy Jey teowa sO, anu 
2.2.3.7 by a fourth-order Runge-Kutta method to yield the 
jet trajectory, width, velocity and temperature decay. 

Because this model is developed only for the 


zone of established flow, the authors try to estimate the 
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length of the flow establishment. This can be done by 
plotting the value of temperature versus distance on log- 
log graph paper and finding the maximum distance up to 
which temperature equals the inlet temperature. Another 
plouren the log-log paper shows thatthe length of the 
flow establishment has a linear relationship with the 
velocity ratio. Another method to estimate the length of 
the flow establishment is by determining the jet half- 
width at the end of the zone of flow establishment. This 
can be achieved by relating the heat flux through the 

end section of the flow establishment to the inlet section 
assuming that the heat loss is negligible in the near 


field region. This gives: 
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or bar= 126° 5b C222 eae) 


where bs is the half-width at the end of the zone of flow 
establishment. 

The value of O45 i also been determined 
graphically by the authors. The value of Cog obtained 


from laboratory and field data is approximately equal to 
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0.5. However, the value of O4 scatters so much that ‘no 
conclusion can be drawn. Besides, no correlation between 
the a, and the velocity ratio is found. 

Paddock, et al (1973) have compared the theo- 
rectical prediction by Motz and Benedict's model to the 
field data obtained at the Point Beach Nuclear Power Plant 
and Palisades Nuclear Power Plant on Lake Michigan, U.S.A. 
They briefly summerize the inadequancies of this model 


as follows: 


1. Buoyant forces, as well as vertical entrainment, 
are ignored in the model development. For den- 
Simetric Froude number less than about 2, buoyant 
Spreading is quite important; for densimetric 
Froude number exceeding about 5, a substantial 
amount of vertical entrainment may be expected. 
Consequently, the types of discharges for wien 
the model may be used is limited. The entrainment 
coefficient, determined from data fitting, must 
actually account for spreading due to buoyancy as 
well as from jet entrainment. 

2. Data presented by Motz and Benedict show consider- 
able variation in the value of entrainment coeffi- 


cient from one situation to another. A particu- 


larly interesting contrast is observed when one 
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compares the Motz-Benedict data for 90° discharge 
to the Romberg-Ayers lake data. Although most of 
the initial conditions are similar, the entrain- 
ment coefficients are different by 10. Limited 
field results indicate that the entrainment co- 
efficient remains essentially constant at a par- 
ticular location with changing values of velocity 
ratio. This behavior is consistent with laboratory 
findings showing a 1 to be relatively independent 
On ues 

The model Ac actually restricted to small ambient 
erent ct The simulation assumes that the jet 
velocity approaches zero at large distances normal 
to the jet centerline. The assumption of Gaussian 
profiles for temperature and velocity and the 
assumption of equal rates of entrainment on offshore 
and leeside of the jet are not valid for ambient 
currents that are not very small. 

The model does not simulate unequal rates of spread 
for momentum and heat. 

The model does not treat turbulent intensity in the 
crossflow, due to the assumption that such turbu- 
lence and its effect on mixing are negligible. 


Although some account of this effect may be made 
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such treatment would be only qualitative in 


nature. 
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Stefan and Vaidyaraman's Model 


Stefan and Vaidyaraman develope a dynamic three- 


dimensional jet model under the effects of wind and 


cross-—-current. 


volume flux 


nb mh 
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ale xe) . - 
heat flux 
nb mh 
H(x, ) Ea (ree Fon Uae oa 
-nb e) P 
momentum flux 
nb mh 5 
M(x_) = GPL “obal eae! 
1? haa te, ‘ Giri aa 


The governing equations are defined as: 


(9624.1) 
dy 
at 


(2.2.4.5 2) 


(252 34.3) 


Parameters m and n in the above integrals are 


finite numbers. 


The authors propose that a value of 
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n=m = 3 is adequate for Gaussian type distributions 
in velocity and temperature. As done in other re- 
search work, the temperature and velocity profiles are 
considered to be similar along the jet axis and are 


described by: 


8 Ae U) Sin e ' 2 , 2 
vu sine Sew oul Ya 7D) 2). @xp (—a (27 7h)e 
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(2.2.4.4) 
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T a Gf 1 h al V 
O a 
C2) 
hy » AL in Equation 2.2.4.5 are known as the simi- 


larity parameters and the authors consider that a value 
of 1.05 is adequate. The assumption of similarity pro- 
files and such a low value of Ay indicate that the 
model results should be applicable for small values 
of Richardson numbers. The authors then apply the en- 


trainment principle with variable values of entrainment 


coefficients to yield the volumetric flux increases in 


the flow direction. 
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The entrainment velocity is assumed to be proportional 
to the difference between the centerline velocity and 
the component of the ambient velocity in the same 
direction. For horizontal buoyant surface jet, the 
entrainment is a function not only of the centerline 
velocity but ee of the degree of stratification, hence, 
a variable entrainment coefficient along the jet must 

be used. The authors adopt an average entrainment co- 
efficient, which is calculated by the following rela- 


tionship, for each cross section: 


Beha ak Apalah (apie atau tsb) (2.2.4.7) 
The horizontal entrainment coefficient dy is 
assigned an average value of 0.059 as in the plane jet. 
The vertical entrainment, ry , through the bottom part 
of the surface jet is inhibited by buoyancy and must be 
equal to or less than horizontal entrainment, depending 
on the overall Richardson number. The authors establish 


the following relationship: 
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The above relationship fits experimental data 
of Ellison and Turner, but only apply with the range 


OnllOaol BER. a hOnO2Siatand GetlGr / A.) < 1.0. 
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The rate of the plume spread is broken into 
two parts, one is due to turbulent mixing and the other 
is due to buoyancy. The authors assume that these 
effects can be superimposed. The horizontal and ver- 
tical spreading due to turbulent mixing are expressed 
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respectively. The lateral spreading of the jet is 
assumed to be the same as that for non-buoyant jet 
except that the absolute spreading angle is distorted 
by the ambient velocity. or is assigned a value of 
0.081 which is considerably lower than the average of 
0.097 taken for non-buoyant jets. Equation 2.2.4.10 
implies that the vertical and horizontal spreading 
ratio is proportional to the square root of the verti- 
Galeand Norizontal mixing cati1o. | 
The rate of lateral spreading due to buoyancy 


is given by the authors as: 
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A value of 0.4 is recommended for C.- As the width of 
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the plume increases, the depth should decrease. Thus: 


As mentioned, the authors consider the effects 
of turbulence and buoyancy on spreading may be additive, 
therefore, the total spread of the jet in the horizon- 


tal and vertical direction are: 
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a ax 
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respectively. 


The rate of change of heat flux is balanced by 
the heat loss to the atmosphere through the water sur- 
face. This gives: 
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in which K is the coefficient of surface heat transfer 


and TE is the equilibrium temperature. 
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Only the wind stress on the surface and the 
drag forces created by the ambient flow are given con- 


Sideration in the momentum flux. The wind force is: 
Pee en Tt (2206 4...16) 
WwW Ww 


The pressure drag exerted on the jet is: 


2 
oma Fi e(U, cos 6) 
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and the friction drag is given by: 
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The drag coefficient Cs Ln seu GLO 23204. 27 
has been given a value of unity. The rate of change of 


momentum flux in terms of the external forces is: 
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in which g is the angle between the wind direction and 
the direction of the ambient current. 

An examination of these equations reveals 
that an important forces created due to density differ- 
ential between jet fluid and ambient has been ignored, 
thus limiting the model's applicability for small 
Richardson numbers only. Besides, the consideration of 
a friction drag when a form drag due to separation has 


been considered does not appear to be realistic. 
2.3 Stolzenbach and Harleman's Model 


Stolzenbach and Harleman use a different 
approach to develop a three-dimensional mathematical 
model for the horizontal surface discharge from a rec- 
tangular open channel. The effect of ambient current 
is also considered. Like in nonbuoyant jet, the authors 
assume an initial core region, void of shear followed by 
the main turbulent region; but unlike most jet analysis, 
the authors divide the structure of each cross section 
into four regions: an unsheared core region, a verti- 
cally sheared region, a horizontally sheared region, and 
a region sheared in both directions, as shown in Figure 
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To derive the basic model equations, the 
authors choose the integral method as in most jet 
analysis. The continuity equations, momentum equations 
in x and y direction, are written for each region along 
with the conservation of heat equation and the trajectory 
equation. The boundary conditions are different for 
different regions. The velocities and temperatures 
at each cross section are presumed to be related to the 
centerline values by similarity functions. The velocity 


and excess temperature are represented mathematically 
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The form of similarity functions are taken 


directly from Abramovich's work: 
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The authors also relate the entrainment 
velocities to the centerline velocity by SBOE SBELaRS 
coefficients. The lateral entrainment coefficient are 
determined by nonbuoyant jet theory alone and the authors 
assume it remains unchanged by buoyancy effects. The 
vertical entrainment coefficient is related to the local 
temperature gradient by the local Richardson number 
using the experimental results of Ellison and Turner. 

The vertical entrainment velocity (We) and lateral 


entrainment (V.) are thus defined as: 


for Z4 = -r-h 

W, = ie - UL cos 8 am OM<s yeov= as 

W o=W Ee) = Ue scosno ah s)< ye< bb) fos G2553816)) 
e e y 1 dx, 

We = 0 Dit Say. 


y ea Wi Ai ‘ Ph if hy . i Be "hi 

a hs SY ; ae , it 7 a 

o i 7 : y a 7 a Aa me 

vs : 7 So ae “in n 


our saa ski eoe | i 
etonivn edt Bie angie. Yaoert# Erte 
ot ,2goeRte yomayoidd yet bocaatond 
‘spol oft of Betelex ab ne istiect rr 
f° seedinen poutasephst feo0i sett es rt 
ee eT anoraitt Ber’ Ww adie * a3 hom 
Leieded bite (a) yeinaten. renee 
ten benitab eu ants we 


®. 


Wie. San RR 9 2 < < 
= cos =i Z n 
e e dx, i 
ay 
A VA £ ( a ) + U. cos © oa algae BuP bc i 2a 
G2 S347.) 
Wie = 0 -r-h > A 


With these additional equations, the authors 
are able to integrate the basic equations over the jet 
cross section to yield a system of simultaneous, first- 
order, nonlinear ordinary differential equations which 
can then Baseaieen numerically. The complete solution 


may be written in non-dimensional form as: 
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In this model, the adequacy of dividing the 
WEG Structure Into four wectanguiar regions ts) an 
question. The presumed jet cross-section structure is 


based upon nonbuoyant jet theory but may not be valid 
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for buoyant jets. The division of the cross section 
also complicates the problem and makes the necessity of 
specifying the inter-regional velocities to fulfill the 
continuity requirement. The form of these interregional 
velocities are unknown and therefore must be guessed. 
The authors also require that no turbulent momentum 
transfer occurs between regions of the jet or between 
the jet and ambient water. This is tantamount to drop- 
ping the Reynold's stress terms in the equations of 
motion. As a consequence, turbulent jet diffusion has 
to be artificially simulated through the entrainment 
coefficient and similarity forms for temperature and 
velocity. This model is also only valid for small 
ambient currents in relation to the initial discharge 
velocity. For significant crossflow, the assumption of 
constant lateral entrainment coefficient may not be 
adequated. The negligence of the pressure drag in the 
governing equations also makes this model valid only for 
small cross-current. If the cross flow velocity is 
sufficiently large, the drag coefficient will no longer 
be negligible. 

Paddock et al (1973) have compared their field 
data collected at Point Beach Nuclear Power Plant with 


Stolzenbach and Harleman's prediction, and have found 
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that the model predicts more rapid decay of centerline 
temperature and velocity. The lateral spread of the jet 
is much greater than the measurements. The model also 
tends to underpredict lateral spread in the region of 
flow establishment. Besides the above comments mostly 
associated with analytical modeling, there are practical 
difficulties in actually obtaining a numerical solution 
to the set of ordinary differential equations. For 

many cases, Paddock et al find that the matrix which 
must be reduced is nearly singular and much precision is 


lost in solving for the derivatives. 


2.4 Phenomenological Model 


2.4.1 Paddock, Policastro, Frigo's Model 


Paddock et al present a phenomenological model 
for the prediction of jet trajectory, centerline temper- 
ature and velocity decay, half-width, half-depth and 
isotherm areas. The model is developed based on the 
measurements made at the Point Beach Unit I and 
Palisades Nuclear Power Plants on Lake Michigan. Velo- 
city and temperature profiles are assumed to be similar 


and of the Gaussian type in the lateral direction. The 
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widths of these profiles are assumed to grow linearly 
with the distance along the jet centerline. After pre- 
liminary examination of the Point Beach data, the authors 
give the following equations: 


For the temperature excess: 


T= TC) (Yq) ; (QeSi.D) 
where 
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in which Ci Cos oy Cy are free parameters. 
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The velocity excess takes a similar form as 


the temperature excess: 
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The temperature centerline trajectory is given in its 


parametric form: 
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The velocity centerline has a form analogous 


to that for the temperature centerline: 
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The free parameters Cy Coy C3, Cys Cs, Ce, 
Co, a Rn Ro K, in the above expressions have to be 
determined by fitting to the field data. With the above 
set of equations, the authors use the minimization 
technique to yield the peste elona. The comparison of 


data is presented graphically. The Point Beach jet 


aie aN eae : a 
K >) a [ a iW i by .) 
i) ne Ms hi 


hy 2 eh ¥ =p ai i 
a 


ist de, Mut | Ae h ‘ie a pai 


it 
r | ne fe a y on i ve so Hh 
oot oi neviw a Bek pails seh a : 7 Ag 
i as ir) mh) i: % re . “ 


ie Ss 


ay 


i '. a 
4 h 7 ae i nes 


Bue oo Cerne aint a pant dni’ Papheaan: eames ont eg 


(vt .8,8 ( @ fr he - ft - (28 
| ° of om 


- ; ¢ | x 
8 hah S) @ a) ao 5 anne 4 eo. ia 
OEsb Sh) a Oe 
Lae apo +59 1g2 sg9 4 (P Presonwrsa 2078 me any 
ts po a 
| od ei eved veeeaneee ony fe a # 7 * rr | 
coved od? d9iw  .cdeb BESkR) Sond wast 4 fal ban 


em an 
“ odd ir belo 


nod seximbai e ae 


rele 


“A. S i 


studies indicate that the centerline excess temperature 
and excess velocity ratios decay at about the same rate 


which can be expressed as: 


1m S Sm x Baer a ji 6 (2.4.1.9) 
ae U6 ak 


Data fitting reveals that generally the 
temperature spreading is wider than the velocity 
spreading in the lateral directions. Comparing their 
data with other established models, the authors comment 
that the Pritchard model compares rather well with 
their limited data but is often conservative when model- 
data discrepancies exist. The Stolzenbach-Harleman and 
Prych models predict too rapid a temperature and velocity 
decay accompanying too great a lateral spread. The 
Motz-Benedict model is too sensitive to the entrainment 
coefficient with little consistent data available for 
its determination for accurate prediction. Because this 
model is developed based on the data taken from Lake 
Michigan, its applicability to other sites still remains 
to be tested. 

Diet ond Elliott and Harkness Model 


Using set of rather extensive measurements 
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in Lake Ontario, Elliott and Harkness choose a pre- 
dominantly phenomenological approach to develope their 
mathematical model for the prediction of the.size, 
shape and orientation of the thermal plume. 19 surveys 
out of the original 53, which were all taken in Lake 
Ontario, are rejected due to instrument failure or 
unfavorable measurement conditions. 

The authors use the least-square analysis to 
determine the functional relationship between the tem- 


perature and the area of an isotherm of the form: 


oer. A 2 (Qa ih) 


Gs eee a ee 
ney Ne 


The term Q(T, - tT) proportional to the heat 
outflow rate and is included to make the expression 
valid for arbitrary values. But its choice of use as a 
scale does not seem to be adequate because the area 
enclosed by the isotherm also depends upon the initial 
momentum, aspect ratio, and the source Richardson Number. 
Before fitting the data to evaluate the regression co- 
efficients Cy and Coe the data are segregated and group- 
ed into three seasons: Winter, spring and fall, and 


summer. It was found that eT LSsequakl to~0Ts0777 073554, 
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0.336; and Cc. has the values 0.712, 0.478, 0.364 for the 
three seasons respectively. The autha’s notice that the 
measured isotherms exhibit a perplexing irregularity in 
shape that would seem to defy any systematization attempts. 
Measurements taken under nominally identical conditions 
produce radically different shapes. The plumes do, how- 
ever, occupy a certain region of the lake surface near 
the outfall and they are elongated in the direction of 
flow away from the outfall. It is, therefore, postulated 
that a Boe factor and not the shape is contained in 
the data. The authors thus define a mean width as: 

2 iss (Dida 23) 


(oo oes ee 
3 12 
(1 + EF.) 


where ¢C..= 0.717 and Cy = 0.552 are the values determined 


S) 
Dy rivrting to. the data. 

Lt is also, difficult ‘to obtain ah expression 
for the trajectory by using the statistical methods 
because the path could be expected to depend on the shape 
of the outfall structure, the outfall velocity, the 
Richardson number, the wind speed and direction, the lake 


current and lake bottom topography. The authors choose 


an analytical method which is similar to the one used by 


Otay a, 


ae 
7 Y 


Fe ‘ee : 


cal wai 


Rw 
A 


wets a poe a 
aig. erty enidos 
niin 2 


4Qmet te suolldiilbedl 


anoiti baad tes ah : iain ‘enell s 
—~<wert oh comply! oar’ vaseedte ae . al fs 
Teer yesh Ne sol ot co) 46 ‘popes eee: 
+o polvoerhb ‘ead. nk Mallnnaats | 
bots fadteoe leseheseah wb a2 ata 
nt. hanipinoo ai edits Nea 


= 
bas 


~. 


($.$.8.8) 


| | | | so a J. i he - 
benisietebh seal ev ogs exe 262,90 = ,0 hae vay ve er 
a at aaa 

AoLeeO gee, M6 MRD OF ahyptis ib AEE, ill 
ebattont sentoinaes hia ahtay ai 


vqera ott ne basaes otf 


| ieee 


| edd ve loptey. w: wel PP) porxatoe a mrad 

Ash odd pans me fini eth, ime ns 

Reams aos she | MAS OPN, ME Be, ON. of 
bina a wai ¥ > oe 

yo Bee eno G98 4 et dokdw Bogiem t MOLINS 


63 


Platten and Keffer (1968). The conservation of momentum 


in the direction perpendicular to ambient current is: 


o— (b ue Sin @) = 0 
i 


(2a too) 


Momentum increased due to entrainment is given as: 


iz ; 1 oe 
oe (OmUl escOS) O)) =" aa. or 
ax 
i . 1 


I, = st} exp (-2(y,/b)*) dy /b) = 972 


Af -1 


Conservation of mass flux gives: 


i= we exp (-(y Vn een) = ¥q_ 
5 L 


In the above expressions, 


to lake current, wind shear, buoyant force, 


(2.4.2.4) 


(2.4.2. 4A) 
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(2.4.2.5A) 


the drag force due 


are all 


assumed to be negligible. Solving *Bquation2y4 32% 3; 
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Equation 2.4.2.6, along with the assumed 
form for the entrainment velocity, can be solved numer- 
ically to yield the jet trajectory. The entrainment co- 
SLiVevenls id Sand) 02 ‘are found by fitting the; data 
to the theorectical solution and turn out to be equal 
to 0.14 and 0.6 respectively. 

Dees Asbury and Frigo's Model 

Asbury and Frigo's study is a completely 
phenomenological approach in attempting to find a re- 
lationship between the plume surface area and the excess 
temperature based upon the field data available from 
other investigators. Seven sets of data, which fulfill 


the criterion set by the authors, are selected in this 
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study. The criteria are: 


a. sufficient data to permit the drawing of at least 
three closed isotherms. 

b. measurements oe the intake and outfall temperatures. 

Cc. measurement of the volumetric discharge flow rate. 


d. measurement of the lake temperature. 


For comparison of the data, the authors plot it in 
a dimensionless form, as suggested by Edinger and Polk 
in their study of thermal plume dispersion, which is 
shown in Figure 2.4. The authors explain that the dis- 
charge Q is chosen as a scaling parameter because it is 
difficult to assign plume-wide values to other variables 
and this prevents the parameterization of scaling area 
in terms of other variables such as ambient velocity, 
ambient diffusivity, and plume depth. The most straight- 
forward parameterization of the scaling area, therefore, 
is of the form is = Q and c is found to fit the data 
best when assign a value equal to one. 

A study on the figure indicates the curve 

fits the data reasonably well. The authors feel that 


the sources of scatter among the data points are: 


a. Large plume-to-plume variation in ambient diffusivity, 
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ambient velocity and plume depth. Although the 
magnitude of the scaling area should be strongly 
dependent upon them, these variables have not been 
included in the parameterization. 

Inaccurate plume areas, especially in the far- 
field region, where plume meandering can produce 
Significant errors in the mapping of thermal plumes. | 
Different outfall geometries. 

The parameterization A ,/2 may not be the best 
choice. There may exist a more suitable para- 
meterization, which would lead to a better 


consistent grouping of the data. 
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CHAPTER .3 


EXPERIMENTAL ARRANGEMENT AND INSTRUMENTATION 


3.1 Experimental Objectives 


The experiments were carried out at the T. Blench 
Hydraulic Laboratory at the University of Alberta. This 
is an experimental study on the heated water discharges 
into a lake with weak and moderate cross-currents. Six 
experiments with different Richardson numbers and cross- 
current intensities were studied. Besides the measure- 
ment of the surface temperature , the temperature field 
in the vertical plane along the jet axis was also recorded. 
The objectives of this study are: 

(jt) ecOnInVeEStigace: thevetrect of cross flow on the 
axis of the jet. 
(ii) to determine the centerline temperature decay. 
(iii) to determine the similarity of temperature 
profiles in both the vertical and lateral direction. 
(iv) to study whether the present theories on surface 
buoyant jets without cross flow can apply to the 


same situation with cross—currents. 


3.2 Experimental Equipment 


3.2.1. Model Basin and Hot Water Supply System 


The model basin ts 50° -— 00" long, 15° - 00” wide, 
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and 2' - 00" deep with a horizontal floor. MThe side walls 
were constructed with 16" x 6" x 8" hollow concrete blocks. 
The inside of the walls was painted to prevent water leak- 
age. A general view on the basin is shown in Figure 3.1. 
At one end of the basin, two identical inlet channels 

were constructed, 4.7 feet apart from center to center. 

One channel was built 2.8 feet from one wall to its center 
and the other was centrally located with respect to the 
basement. These inlet channels are 5 feet long, 1.4 feet 
wide and 2.0 feet deep. The thermal effluent was supplied 
by a water heater which has a burner input of 900,000 BTU 
and coupled with a pump, has a maximum capacity of 60 gpm 
at a constant temperature of up to 130°F. (Figure 3.2 shows 
the water heater). The heated water was discharged into 
the inlet chamber from the bottom and then passed over a 
V-notch to the inlet canal before entering the basin. The 
V-notch was used to measure the flow rate. A 2" diameter 
telescopic pipe for overflow was fixed in the inlet chamber 
to maintain a constant head over the V-notch. The invert 
of the inlet canal was 14" from the floor. The schematic 


diagram of the whole system is shown in Figure 3.3. 


3.2.2 Cross-Flow 


The cross flow was created by a series of closely 
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spaced rectangular jets. This jet effect was created by 
cutting seventy-nine rectangular slots of 4" x 1/16" ona 
5" diameter plastic pipe which was 20 feet long. The 
spacing between slots was 3" from center to center. Those 
slots were completely submerged with the top coinciding 
with the water surface. The axis of the pipe was about 

4" away from the inside of the wall and the jets were 
expected to merge together at a distance of about 10 inches 
from the points of discharge. The pipe was connected to a 
pump which was submerged in the sump. The capacity of the 
pump was 0.1 cubic feet per second. Figure 3.4 shows the 
cross flow system. In some of the later runs, the pipe 
from the sump to the cross flow pipe was wrapped with an 
electric heater wire to increase the temperature of the 


cross flow. 


3.2.3 Outlet Water Collecting System 


The outlet water was collected by two separate 
systems. At the end of the basin, the water flowed over 
3' - 7" long and 8" high tailgates into the outlet chamber. 
The tailgates were constructed from marine plywood and 
rested on 8" high and 6" wide concrete blocks fixed on 
the floow. A 6" diameter pipe was fixed at the bottom of 
the concrete block to drain the basin. Another pipe out- 
let from the outlet chamber led the water collected to the 


sump. 
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me 


The other collecting system (Figure 3.5) consis- 
ted of 4 sections of 5 feet long plastic pipes. Nineteen 
slots 4" x 1/16" were cut on each section and the spacing 
between slots was 3" from center to center. These pipes 
were fixed at such a height that the slots were completely 
submerged. These four sections of pipes were jointed by 
one brass pipe which extended to the outlet chamber, but 
each section had its own valve to control the amount of 
pucevarn Those pipes were fixed on the side of the basin 
opposite to the pipe that generated the cross-current. The 
purpose of this system was to collect the extra discharge 


from the cross-current to maintain a cross-flow. 
3.3 Instrumentation for Thermal Measurements 
3.3.) Thermistor Probes 


Thermistor probes were used to measure the tem- 
perature of the heated water in the basin. Thermistors are 
resistors with a high negative temperature coefficient of 
resistance. As the temperature increases, the resistance 
goes down and vice versa. They are semi-conductors of 
ceramic material made by sintering mixtures of metallic 
oxides such as manganese, nickel, cobalt, copper, iron 
and uramium. Their electrical characteristics can be con- 


trolled by varying the type of oxide used and the physical 
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size and configuration of the thermistor. 

The thermistors used in these measurements were 
glass-covered probes ( code number GB 32P2) supplied by 
Fenwal Electronics, Framingham, Massachusetts. They are 
2 inches in length, 0.012 inches in diameter and the lead 
material is Dumet. They can be used satisfactorily up 
to 1000°F, but the manufacturer suggests that for maximum 
stability, operating temperature should not exceed 600°F 


which well covered the temperature range in our experiments. 
Deed <i Probe Calibration 


The calibration of the probes was done in a water 
bath. Thermistors were submerged in water and the water 
temperature was measured by a quartz-thermometer (model 
HP-2801A). Such a thermometer monitors the frequency 
Changes of quartz crystal oscillation. Because the quartz 
thermometer converts temperature into frequency rather than 
into resistance or voltage, it is free of the problem of 
lead resistance and noise pickup inherent in other types 
of temperature measuring systems. The temperature of the 
water was varied and the response of the thermistors was 
read from the digital voltmeter (model HP-3480B). Non- 
linear least square technique was adopted to fit an equa- 


tion, which converted the voltage reading into temperature. 
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Atypical calibration curve, for the themistone is given in 
Figure 3.6. The calibration equations are given in 
appendix. 

1S es, Location of themistor probes 

Eighteen probes were positioned at different 
locations to measure the temperature. Sixteen of them 
were mounted on a rake trolley (Figure 3.7) but the first 
and the last probes were not connected to the data acgi- 
Sition system. The rake trolley was powered by an electric 
motor to move in the vertical and lateral directions. 
Safety limit switches were installed in both directions to 
avoid any possible damage to the system. The glass probes 
were mounted on brass tubes, 1.5 feet long, spaced 3 inches 
from center to center. Two 1.6 feet long solid brass tubes 
were also mounted on the same rake to avoid damage to the 
probes due to accidental impact on the bottom of the basin. 
In order to locate the probe system along the jet, the rake 
trolley was mounted on another trolley that could be moved 
manually in the longitudinal direction. One probe was 
mounted at the inlet and the other was positioned by the 
side wall, and they were used to measure the inlet and 
cross flow temperatures respectively. The location of the 
rake was sensed electronically and the coordinates were 
punched on the paper tape along with thermistor probe 


readings. 
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3.3.4 Data Acquisition System 


The Data Acquisition System consisted of the 
following elements: | 
(1) Co-ordinate positioner 
(2) 8000 A Digital Multimeter (John Fluke MFG. Co) 


(3) 100 channel scanner - Digitec model 635 
(United Systems Corporation) 


(4) Scan counter 
(5) HP-3480 B Integrating Digital Voltmeter 
(6) HP-3482 A D.C. Range unit 


(7) Model 623 Punch Controller 
(United Systems Corporation) 


(8) Papertape punch (Franklin electronics, inc.) 


Figure 3.8 shows a general view of the system. 

Each thermistor was assigned to a specified channel 

scanner in addition to the X, Y and Z co-ordinates. Change 
in temperature would result in the change of the resistance 
of the thermistor probes and eventually change the voltages 
which were fed into the channel scanner. The scanner then 
routed the analogue signals to a converter where the sig- 
nals were digitized and converted to binary coded decimal 
(BCD) form. The punch controller received the BCD signals 
and converted them into ASC1l code before routing it to the 


papertape punch. The papertape was read by the 1145 FEBD 
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processor tape reader where the ASC11 code was automatically 
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translated to EBCDIC code, the code used by the IBM 360 
computer, and stored in the computer files. A pre- 
written program was then used to calculate the temperature 
from the input file and printed out along with their aver- 
age value. Figure 3.9 is a flow chart of the Data Acqui- 


Sition and processing systems. 
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CHAPTER 4 


EXPERIMENTAL RESULTS AND ANALYSIS 
4.1 General 
In this study, six experiments with Richardson 


HuUMDerowOL Ol0S57 OO, lr +0..29) 0256," 0l7 > rand” bal! were 


carried out. The source Richardson number is defined as 


the inverse of the square of the source densimetric Froude 


number : 
1 ay 
R. — 2 = ( os 8 Be) FPS) gh, 
10 (Fo) Dipped ST en 
UG 
where he jet depth at outfall 
OS water density at cross flow temperature 
0 water density at the outfall 
O 
us jet velocity 


A previous study by Pande indicated that the 
behavior of surface discharge is governed strongly by the 
Richardson number. Though the exact values that divide 
the Richardson numbers into groups of large, moderate or 
small R;, are not well defined, 0.1 and 0.7 are believed 


to separate low Richardson number from the moderate, and 
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the moderate from the large Richardson numbers respectively. 
The rate of expansion of the surface discharge in any one 
Richardson number group is very different from that in any 
otter group. Pande found that the rates of spreading and 
decay of excess temperature for surface discharges without 


cross flow as follows: 


a. Small Richardson number 
1 al 
U. o = ne = bex hex 


b. large Richardson number (assuming re) 


U ! ay : b 
o ax ore ee f Ox 
m a1 /A3 m eV 
c. moderate Richardson number (assuming hex“) 
1 1 4/3 S 
C= 0 Ura. eis Ty * 2/3 bex hex 
1 i 1d 
qee=i1 Uz Stas Tc B57 3 Dex h “x 
ee My 3s 
where Ue centerline velocity 
US kate centerline excess temperature 
b lateral distance from centerline where 
AD = AD 
m 
D4 longitudinal distance measured from inlet 
h distance measured downward from surface 


where Ain = AT 
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This study presents an investigation on the 
effect of cross-current on the surface jet behavior. The 
Reynold number was varied from 2840 to 14100 and the sur- 
face discharges were turbulent in nature. The details of 
initial conditions for all experiments are given in 


table 4.1. 


4.2 Temperature Field 


The temperature was measured by the thermistor 
probes as described in chapter 3. Experiments with 
Richardson numbers of 0.75, 0.1, 0.035 were run using the 
inlet 2.8' away from the side wall. The inlet had a bah. 
ratio of 1.19. In these three runs, the measurements 
started at a distance 0.5' from the inlet and the subse- 
quent sections were 0.5 feet apart from each other. Other 
experiments with Richardson numbers of 0.29, 0.58, 1.14 
were run using the centrally located inlet which had a 
Do/hs racic of 1.18. The subsequent cross-sections were 
spaced 0.25 feet apart. In each experiment, a preliminary 
analysis on the surface temperature distribution was car- 
ried out to determine the jet axis, and then the tempera- 
ture! profiles’ onthe vertical plane along the axis were 
measured at an interval of one foot. The reliability of 
the temperature readings was checked by calibrating the 


thermistor probes after the experiments and the agreement 
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Was found within 0.2 degree centigrade of the initial 


readings. 


43 Cross-flow Velocity BPeor lee 


The velocity structure of the cross-flow 
(Figure 4.1, 4.2) was measured by Mr. Brent Berry using the 
Hydrogen-Bubble technique. Because the velocity was not 
uniformly distributed in all regions, average value was 
used when relating the parameters within that particular 


region. 


4.4 Excess Temperature Profiles 


The excess tempeyature i Tete) profiles 
measured on the surface yield a jet-like behavior with 
maximum excess temperature at the axis of the jet and de- 
creasing towards the edges. The profiles also indicate 
clearly that the excess temperature profiles are not sym- 
metrical about its axis due to the pressure and shear drags 
created by the cross-current. A typical set of cross- 
section profiles is shown in Figure 4.3 for R;, = 0.035. 
The excess temperature profiles for cross-sections perpen- 
dicular to the jet axis for the same Richardson number are 
shown in Figure 4.4. Measurements along the axis in the 
vertical plane show the maximum excess temperature is at 


the surface and AT decreases towards the bottom. These 
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profiles for all cross-sections are similar along the jet 
axis. Typical set of vertical excess temperature profiles 
is shown in Figure 4.5 for Rio 0.75. The presence of 
the cross flow does not seem to affect the shape of the 
profiles. The excess temperature profiles in the vertical 
direction of all runs show the same pattern regardless the 


Richardson number. 


Gas Similarity of Surface Temperature Profiles 


For all runs, the temperature excess profiles 
Bene checked “fon similarity by plotbingg(Te— 0 97/{t ) =) 
versus Y,/b,, for the upstream and downstream sides 
separately. The profiles have been found to be similar 
and described by the Gaussian Curve. A typical profile 
is shown in Figure 4.6 for Riana Las 

The exponential (or Gaussian) curve 


as exp (- 0.693(-1)%), and the Schlichting-type distri- 


AT bm 
2 


5 Be a 
ae AT Set 80 7 (2) 2) are plotted on the same 


AT be 


figure for comparisiton™ Both the exponential curve and 
Schlichting-type distribution agree well with the data up 
to n = 1.0; but the exponential curve gives values slightly 
below those obtained from measurements after n= 1.0 eee 


seems the best fit to the data is the Schlichting-type 
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distribution but the exponential curve has the advantage 
of being simple to work with in integral operations. The 


existence of the cross flow does not affect the similarity 


of the excess temperature distribution on the surface. The 


Similarity plots for other runs are given in Figure 4.7 to 
Figure 4.11. 

Though the assumption of similarity still held 
for each experiment, the shape of the similarity profiles 
on the upstream is affected strongly by the cross flow, 
the exponential curve being obviously a very poor approxi- 
mation in the chosed neighbourhood of n#=1.0. The down- 
stream profiles are not affected by the cross-current at 


all and can be treated as in the case of zero cross flow. 


4.6 Similarity of Temperature Profiles on the Vertical 
Plane along Axis of Jet 


The excess temperature profiles on the vertical 
plane were also found to be similar for all sections along 
the axis of the jet for all runs. A typical non-dimen- 
sional plot of AT/ AT versus Z,/h fore Ra) 10-0255 1s 
shown in Figure 4.12 

The Schlichting type distribution fits the data 
best for Z,/n< 1.0 but lies above all experimental data 
for Zi/h > 1.0. The exponential curve describes the data 


reasonably well. In the case of no cross-current, Pande 
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found the Schlichting type distribution, which predicts 
higher values than the exponential curve towards the bottom, 
agrees well with his data. Comparison of these two cases 
indicates that the excess temperature decays much faster 
towards the bottom in cross flow than in the case of zero 
cross flow. The similarity plots for other Richardson 


number are given in Figure 4.13 and Figure 4.14. 


4.7 Excess Temperature Decay along ythe Axis 


The decay of the excess temperature for all runs 
is plotted in Figure 4.15 and 4.16. The decay of the ex- 
cess temperature is not linear but decrease sharply at the 
beginning and then more gradually along the axis. The 
relatively distinct locations of the data indicate the 
influence of other factors such as Richardson number and 
aspect ratio. 

A dimensionless plot of He 7 AT. against 
x4 / YAS. where AT ia is the average of the upstream and 
downstream excess temperature, is shown in Figure 4.17. 
The length scale has been chosen as YAS following the 
work Of Pani/ (1972) on bluff wall jets. 

In Figure 4.17, the data obviously can be divided 
into-two-groups with>-one- located -above-the-other.. The 


upper group contains the data obtained by using the side 


inlet and has a smaller velocity ratio. The lower group 
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of data points are collected using the centrally located 
imtet. and a much higher velocity ratio. sit 1s sditeicult 
to explain the phenomenon that the measurements with a 
stronger rare eur ere do not decay faster than those with 
a weaker current. Also, it appears not very logical that 
the data with stronger current is closer to the data 
obtained by Pande without cross flow. Though Carter's 
laboratory experiments have shown that the center line 
temperature decay has a form apparently independent of 
Ene velocity ratio,, we believe, that thesincrcasing amount 
of entrained fluid due to a stronger current should 
increase the dilution. Figure 4.17 also shows that all 
the data lie below the data of Pande. The faster decay 
rate is due to the existence of the cross-current which 


increases the rate of entrainment and dilution considerably. 


4.8 Correlation of the Temperature Coefficient with the 
Velocity Ratio 


For each run, AT,/ AT 1q 18 Plotted against 
Spt VA, and the equation that fits the data best is 
determined. The equation assumes that AT, is a function 


Cine 2, vA. 


O and ieee and all the other variables are 


included .in. a coefficient oe a is then plotted against 


the velocity ratio 4 = U,/U;- 
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Large Richardson Number 


Previous study by Pande indicates that the excess 
temperature on the axis varies as the -2/3 power of the 


distance along the axis, and can be written as: 


AT na -1/3 al 


= -2/3 
“aT, SS OG, | a + Cy) 


Pande evaluated this equation using his experimental data 


and gave: 


as = x -2/3 
—tS = 3.1664 (R; .) eee psa) Sy) (4.8.1) 
AT, Ao 
B72 
The plots of ( Ae We) versus x,/ VB scone ae Ono 


and 1.14 are given in Figure 4.18. Table 4.2 A gives the 
summary of the calculations. 


The best fit equations determined from the data 


Cine). 

for Rig = 10.7.5 
a s027 ce ee ae ay Yea (4.8.2) 
AT, vie 

for Rio = 1.14 
—ma = 51684 Re ee 4 oe ee (4.8.3) 
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Figure .4.21°A. 1s the plot of C," against: t/a 
for large Richardson numbers. Because of the limited data 
available, no equation can be written. It only shows the 
trend that Ct Pane to a minimum as the velocity fatie 
approaches an optimum value and after that, Cy increases 


with decreasing a. 
Moderate Richardson Number 


For moderate Richardson number, Pande found the 


equation: 


“nna ( ae “ 1 Be ee (4.8.4) 
=A 29s (R —— +1. ae 
AT, Lo VR 


In our study, the equations are found to be: 


TO. Neo = On 2 9 


AT 1735 a) -2/3 
feel ise (R. eee es 1g) (4.8.5) 
AT. Vi, 
for Ri, = 0.56 
ATma 1/3. % -2/3 


ya+Qie, Di) (4.8.6) 


= le3o. (RB: 
AT> ( me) ( a. 


respectively. The summary of the calculation is given in 


/2 


3 
Table 4.2 B and the plots of ( ays AT) versus 


4 tbo for the two Richardson numbers are given in 
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Figure 4.19. From Figure 4.21 B, the temperature coeffi- 
cient C, shows the same trend as in the large Richardson 
number, but the optimum value foro is different. 
Small Richardson Number 
For small Richardson number, Tamai, Weigel and 
Tornberg give the equation 
Ae . a 
og 
a os 7.0 - ( —) (4.8.7) 
for temperature excess decay for a circular jet without 
Cross flow. (D, in Equation 4.857 4s the diameter of jet 
outlet and X is the longitudinal distance). Our experi- 


ments yield the equations 


potas = Ones 
AT x 2 
en 5 eo! er) (4.8.8) 
fore, = 0.035 
10 
AT x 2a 
me = 5.61 (—= + 1.9) (4.8.9) 
AT VE 


It seems that in the small Richardson number 
group, the Richardson number has no effect on determining 


the temperature decay. The two equations determined from 
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our data differ only slightly from each other in the 
temperature coefficient. Table 4.2 C is the summary of 
the calculation. AT. /, ATins. 28 plotted against %/V AQ 
in Figure 4.20 for R;, = 0.1 and 0.035 respectively, and 
on is plotted against w/o fin Figure 4,21) (ce 

A comparison on the plots for different 
Richardson number groups shows that the temperature co- 
efficients for small Richardson number are much greater 
than the other two. groups. The data show the same trend 
as C, decreases to a minimum for an optimum value ® and 
then increases as a is decreasing. The rate of change 
Of @, with respect cotyl/o is different foriditterent 


t 


Richardson number groups and so is the magnitude of C.. 
4.9 Isotherms 


The surface isotherms for all experiments are 
plotted in Figure 4.22. Observation of the isotherms 


snows thet for Ri, = 0.29 and 0.75, the surface jets 


spread more on the downstream than on the upstream. This 


behavior is expected because of the additional pressure 
drag exerted on the upstream by the cross flow. For 
Rig = 0-035, the width on both sides are approximately 


equal; but for runs with Reo” La i4, Of and 025657 we 


find that the jets spread more on the upstream than the 
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Table 4.2 Values of Temperature Coefficient 
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FIGURE 4.21 VARIATION OF TEMPERATURE COEFFICIENTS WITH VELOCITY RATIO 
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FIGURE 4.22 SURFACE ISOTHERM 
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downstream side. Hoopes et al. did not observe such 
behavior in their experiments. Policastro et al. stated 
that large eddies on the down-current side should cause 
more extensive spreading on that side. The exact reason 
for this phenomenon is not known but we believe that when 
the cooler cross flow reaches the front of the surface jet, 
it dives underneath the jet simply because of its heavier 
density. As the ambient current passes under the jet, 

the motion creates a shear force between the two layers 
and eventually sets up a circulation within the jet which 


might cause a greater spreading on the upstream side. 


4.10 Jeti Trajectory 


The jet trajectories for all experiments are 
shown in Figure 4.23 and Figure 4.24. Three regions can 


be observed: 


1. Undeflected region where the axis of the jet 
coincides with the centerline of the inlet. 

2. Region of maximum deflection where the cross flow 
affects the jet to the maximum and deflects the 
axis of the jet towards downstream at a very sharp 
race, 

3. Region where the jet is directed along the already 


deflected direction. 
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FIGURE 4.23 JET TRAJECTORIES ON SURFACE 
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These three regions are shown clearly for all 
runs except for the run with Rio = 0.29. In that parti- 
cular run, the strength of the cross flow is so small as 
compared to the jet velocity that it does not affect the 
axle atfall. In other runs, the degreerof bending of the 
axis depends strongly on the strength of the cross-current. 

The angle of deflection © is defined as the 
angle between the inlet centerline and the axis of the 
deflected jet. Values of tan 6 for all experiments are 
given in Table 4.3 and are also plotted against the velo- 
city ratio «in Figure 4.25 along with the data from Motz 
and Benedict. Though the data scatters quite a lot, a 
general trend is noticed. The same data is again plotted 
on a log paper in Figure 4.26. In this plot, the angle 
of deflection is normalized by the angle of discharge. 

With the additional data from Motz and Benedict, a straight 
line can be drawn through the data points. To observe 

the effect of Richardson number on the angle of deflection, 
we plot tan © (eee against o in Figure 4.27. A straight 


line having the form 


Tan © = (R,Q) (0-741 = 0.045 ) (AS10010 


can be drawn through the data. k is an exponent that has 


to be determined by a trial and error procedure; and is 


found equal to -1/3 for small Richardson number and 3/4 
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for both the moderate and large Richardson numbers. 
45110 Variation of Half-width 


The surface temperature half-width b is defined 
as the distance on the surface measured perpendicular to 
the jet axis where the temperature excess (T - Dep is 
half the axis temperature excess fee = Danie dimensionless 
plot of b/b, vs. 7 De is shown in Bigure 4.26." Jthe data 
points fall within a kidney-shape area except the upstream 
half-width for Rae = 0:56. Tn othat. particular run,. the 
upstream half-width grows very rapidly and has a magnitude 
much greater than the other measurements. The half-width 
for upstream and downstream are plotted separately for all 
experiments in Figure 4.29 and Figure 4.30. For large and 
small Richardson numbers, the half-width grows linearly 
with ae but for moderate Richardson number, it varies as 
Soe This behavior is consistent with previous study by 
Pande without ambient current. At some distance from the 
inlet, the half-width drops off which indicates that the 
ambient current comes into effect therein and restricts 
the spreading of the jet towards the upstream. The up- 
stream half-width of run with R;, = 0.75 shows a very 


irregular pattern and no reasonable explanation is found. 
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Large Richardson Number 


For large Richardson number, the half-width 


varies as x, A non-dimensional plot of b,/b, against 


pi7D6 for Rig = 0.75 and 1.14 is shown if Figure 4.31. 


where be is the average of the upstream and downstream 


half-widths. Pande gave the following equation in his 


study. 
2b x 
eS! WS ae I 


A half-width coefficient Ch is defined as: 


2b 
a 


b 


78 ra 
=C, (Rio) Bo + Co) 


In Figure 4.31, the straight lines have the equations: 


for Ri, = 0.75 


The 


1/30 
= 1.01 (Rio) Ce irae) (4.11.2) 
O 
[amis bee 
1/a x, | 
= 2.09 (R;,) ieee) (4.11.3) 
Bo 


tabulated values for Cy and a are given in Table 4.4 
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and the plot of CG, against 1/a is shown in Figure 4.34. 
Moderate Richardson Number 
For moderate Richardson number, the half-width 


coefficient Ch is defined by the following equation: 


2b 1/6 /3 


ope ep CuO) 


a A 
( 7 + Cc) 
re) O 


Pande found the following equation fit his ex- 


perimental data well: 


2b 1/6 cial 4/3 
bo ~ 0.81 (Rig). iy. * 2) (464054) 


B74 
(2b. /b,) is plotted against x1 /d . invPigurer4.32 for 


ee 0.29 and 0.56. The equations that fit the data are: 


OL Rio = 0.29 


2b, 1/6 a 4/3 
Soret = 1,13 (Re? nr (45a) 
O O 


and ror Ris = 0.56 


2b x 4/3 
ae 136 Qin go 5) (ALG) 
bo O 


The parameter Cy is determined from the experimental data. 
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Small Richardson Number 


For small Richardson number, the half-width is 


duvecte ly proportion tox Pande has no data for small 


ee 
Ri but we define the half-width coefficient as: 


Wiegel and his associates related the half-width 


of the jet to the Richardson number and longitudinal dis- 


tance for a circular surface jet. Their findings give the 
equation: 
2b 1/8 x 
Be lA (Ra ,) / yale (42 1a) 
b b 
oO Oo 


The non-dimensional plots for small Richardson 
number are given in Figure 4.33. The forms of the equa- 
tions we found are: 


for RR. = 0.1 
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The values of Cy, nor all ‘experiments sare tabu— 
Jated in Table 4.5. and then plotted against 1/a in 
Figure 4.34. We find that Cy, increases to a maximum value 
with ® approaching an optimum value and then further de- 
crease in ® results in decreasing values of Cc): The 
optimum value for q and the maximum value for C, are 
different for different Richardson number groups. The 
ratio of the upstream (and downstream) half-width to the 
total half-width against x, is plotted in Figure 4.35. 
The plot shows that the ratio of upstream (and downstream) 
half-width to the total width is a constant regardless the 
distance along the jet axis for all experiments. This 


ratio is plotted against the velocity ratio in Figure 4.36. 


4,12 Variation off Half—Depth 


The experimental results on the half-depth for 
all runs are shown in Figure 4.37 and Figure 4.38. For 
large and moderate Richardson numbers, the half-depths 
are practically invariant with the axial distance as pre- 
dicted by Pande. For small Richardson number, the theory 
predicts a linearly increasing half-depth. In the run with 
Rio = 0.1, we observe a decreasing in half-depth along the 
axis; in other words, the jet is getting thinner and thin- 


ner. The reason for this abnormal behavior is not well 


understood but we believe that the buoyancy plays an 
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Table 4.4 Values of Spreading Ratios 
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Table 4.5 Values of Half-Width Coefficient 
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important role even at this value of Rage eh corresponding 
dimensionless plot is shown in Figure 4.39. 
In order to determine the effect of cross flow 


om he depth scale, we detine a depth coerricient '¢. sac: 


-1/8 
h = 

for large and moderate Richardson number. Table 4.6 gives 
the values of on for different velocity ratio. Figure 
4,40 is a plot of Cy against 1l/fa . Cc, can be taken 
approximately as 0.3 for all 1/a within the range of 


large and moderate Richardson numbers. 
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Table 4.6 Values of Half-Depth Coefficient 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 
FOR FURTHER RESEARCH 


Conclusions 

Thermal discharges from power plants are capable 
of causing considerable damage to our environment. A 
thorough understanding of the behavior of the surface 
discharge is necessary to minimize its impact on the en- 
vironment and also to maintain the efficiency of the power 
plant cooling system. 

Based on the results of six experiments (with 
two experiments each for small, moderate, and large 


Richardson number groups) with cross-currents, the follow- 


ing conclusions are drawn: 


The excess temperature profiles in the lateral 
and vertical directions have been found to be similar and 
could be described either by the exponential curve or the 
Schlichting-type distribution. But the excess temperature 
profiles on the upstream side, even though similar, are 
strongly distorted By the “cross-current; and) neither the 
exponential. curve: nor, Schlichting-type, distribution, could 
describe them satisfactorily. The decay of the centerline 


excess temperature is affected to some extent by the 
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strength of the cross-current. The temperature generally 
drops off very rapidly within a short distance from the 
inlet and then decreases gradually to the ambient tempera- 
ture. The jet trajectories depend strongly on the velocity 
ratio. Though the symmetry of the jet is affected by the 
existence of the cross-flow, the ratio of the jet spread- 
ing on either side of the jet centerline to the total 

width of the jet at any cross-section remains approxi- 
mately invariant in the axial direction. This ratio de- 
pends on the velocity ratio as well as on the Richardson 
numbers. Further, it has been found that: 


for small Richardson number: 


for large Richardson number: 


3 
es rae De xy h « oT 


and for moderate Richardson number: 

4/3 

ee b pa h 5 

AT. « Seyi rr « x) 
qx 


The coefficients in the above equations have 


been found to be functions of the velocity ratio. 
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Suggestions for Further Research 


The surface spreading of the heated jet with 
cross-currents should be studied in more detail with 
measurements of the velocity field. Further, for each 
Richardson number group, a few more velocity ratios should 


be investigated. 
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Calibration Equations 
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Table A-1 Calibration Equations for Thermistor Probes 
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Measured Excess temperature Decay ("C) Along Jet 
Centerline on Surface 


Table A-2 
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